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Protein phosphorylation homoeostasis is tightly controlled and pathological conditions 
are caused by subtle alterations of the cell phosphorylation profile. Altered levels of 
kinase activities have already been associated to specific diseases. Less is known 
about the impact of phosphatases, the enzymes that down-regulate phosphorylation 
by removing the phosphate groups. This is partly due to our poor understanding 
of the phosphatase-substrate network. Much of phosphatase substrate specificity is 
not based on intrinsic enzyme specificity with the catalytic pocket recognizing the 
sequence/structure context of the phosphorylated residue. In addition many phosphatase 
catalytic subunits do not form a stable complex with their substrates. This makes the 
inference and validation of phosphatase substrates a non-trivial task. Here, we present 
a novel approach that builds on the observation that much of phosphatase substrate 
selection is based on the network of physical interactions linking the phosphatase to the 
substrate. We first used affinity proteomics coupled to quantitative mass spectrometry to 
saturate the interactome of eight phosphatases whose down regulations was shown to 
affect the activation of the RAS-PI3K pathway. By integrating information from functional 
siRNA with protein interaction information, we develop a strategy that aims at inferring 
phosphatase physiological substrates. Graph analysis is used to identify protein scaffolds 
that may link the catalytic subunits to their substrates. By this approach we rediscover 
several previously described phosphatase substrate interactions and characterize two new 
protein scaffolds that promote the dephosphorylation of PTPN11 and ERK by DUSP18 and 
DUSP26, respectively. 
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INTRODUCTION 

Protein phosphorylation is a common post-translational modi- 
fication governing signal propagation (Mann and Jensen, 2003). 
The concerted activity of kinases and phosphatases modulate 
protein phosphorylation levels and control key physiological 
processes, such as migration, proliferation, inflammation, and 
apoptosis (Graves and Krebs, 1999; Manning et al., 2002a,b). 
Till recently protein phosphatases have been considered unin- 
teresting housekeeping enzymes and have received less attention 
compared to kinases (Bardelli and Velculescu, 2005). However, 
evidence accumulated over the past decades have indicated that 
this enzyme class plays an important regulatory role and that 
the deregulation of the concentration or activity of specific phos- 
phatases correlate with a variety of human disorders (Wera and 
Hemmings, 1995; Tonks, 2006). Notably, approximately 40% of 
protein phosphatases are implicated in tumor development, high- 
lighting the central role of this enzyme group in growth regulation 



and identifying some members of this enzyme class as promis- 
ing therapeutic targets (Julien et al., 2011; Liberti et al., 2013). 
One of the problems in the characterization, on a large scale, of 
the functional role of members of the phosphatase family is the 
lack of a simple, robust, method to identify physiologically rele- 
vant substrates. Many phosphatases have low intrinsic enzymatic 
specificity and are able to de-phosphorylate many substrates non- 
specifically in vitro (Tremblay, 2009). Alternative methods such 
as the use of trapping mutants (Blanchetot et al., 2005) are often 
used, but the identification of direct phosphatase substrates still 
remains a challenge. 

In order to characterize new modulators of some key cancer 
associated pathways and to identify their direct targets, we have 
recently proposed a novel strategy based on a phosphatase high 
content siRNA screening combined with modeling and simula- 
tion. This approach enabled the identification of 62 phosphatase 
catalytic or regulatory subunits whose down-regulation affects 
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one or more of five readouts linked to cell proliferation: ERK, 
p38, and NFkB activation, rpS6 phosphorylation and autophagy 
(Sacco et al, 2012a). However, this approach was not designed 
to identify the direct phosphatase substrates, responsible for the 
phenotypic effect. 

Here we delineate a strategy to identify protein scaffolds that 
may contribute to substrate recognition specificity by bridging 
the phosphatases to their targets. To develop this strategy we 
focused on eight phosphatase subunits whose down-regulation 
was found to affect ERK and/or RPS6 phosphorylation and are 
therefore modulators of the RAS-PI3K pathway (Figure 1). To 
identify new phosphatase substrates involved in the control of 
the RAS-PI3K pathway we first built a protein interaction net- 
work (PPI) by combining information extracted from protein 
interaction databases and results from new affinity purification 
experiments. This analysis confirmed that the identified phos- 
phatase interactors often act as molecular bridges linking enzymes 
to substrates. In addition we independently validated a subset of 
these predictions 

RESULTS 

THE PHOSPHATASE INTERACTOME 

We have used the results of the siRNA screening (Sacco et al., 
2012a) to select eight phosphatases that modulate the activity of 
the RAS-PI3K pathway. The phosphatase catalytic or regulatory 
subunits were cloned in frame C-terminal to an SF-TAP cassette 
and transiently transfected in HeLa cells. These constructs direct 
the synthesis of four tyrosine phosphatases (PTPN21, PTPN3, 
DUSP18, and DUSP26), three components of the PP2A holoen- 
zyme (the PPP2R3C regulatory subunit, the PPP2R1A scaffold 
subunit and the PPP2CA catalytic subunit) and the PPP3CA 
(calcineurin) serine/threonine phosphatase. As negative control, 
HeLa cells were transiently transfected with the empty vector 
SF-TAP. Since in our siRNA screening phosphatases controlling 
the activity of the RAS-PI3K pathway were identified in HeLa 
cells stimulated with TNFa for lOmin, we decided to perform 
the affinity purification experiments in the same experimental 
condition. Thus, phosphatase transfected cells were stimulated 
with TNFa for lOmin or left untreated. While the control cells 
were grown in a medium containing natural amino acids, phos- 
phatase transfected cells with or without TNFa were, respectively, 
grown in media containing isotopically labeled lysine and argi- 
nine amino acids (SILAC) (Ong et al., 2002). After lysis, phos- 
phatases, and regulatory subunits were affinity purified and ana- 
lyzed by mass spectrometry (Table SI), as described in Materials 
and Methods (Figure 2). Contaminants binding to these baits 
were identified by their equal abundance in both (Table S3), 
the affinity-purified phosphatase sample and the negative con- 
trol (Meixner et al., 2011), whereas true co-purified interactors 
to a given phosphatase were identified by selective enrichment 
of their peptides (Table S2). Only those that were significantly 
enriched in our samples were considered for further analysis, 
as described in Materials and Methods. As shown in Figure 3, 
this strategy resulted in a highly connected interaction network. 
Approximately 10% of the identified interactions have already 
been reported in the literature. Indeed we were able to recapit- 
ulate most of the interactions occurring between the catalytic 



and the scaffold subunits of the PP2 holoenzyme, which, as 
expected, share a significant number of common interactors, 
many of which are regulatory subunits. These observations taken 
together with the validation by coimmunoprecipitation assays of 
some of the newly identified phosphatase interactions (Figure SI) 
confirm the reliability of our experimental approach. In addi- 
tion in Figure SI, we demonstrated that our affinity purification 
experiments enabled the identification of new phosphatase inter- 
actors (the dynein protein, DLC1 and the serine threonine kinase, 
ATM) already involved in the regulation of the autophagy pro- 
cess. With the exception of the PPP2CA-DLC1 binding, both 
DUSP26-ATM and PTPN21-DLC1 associations are decreased 
upon an autophagic stimulus (starvation), suggesting that these 
interactions may have a regulatory role in the autophagy process. 

Since the affinity purification was also carried out with or 
without incubation with TNFa we can also provide dynamically 
regulated interactions in response to TNFa treatment. As shown 
in Figure 3, a few interactions are positively (green edge) or nega- 
tively (red edge) regulated by TNFa incubation. The vast majority 
of the co-purified ligands, however, are TNFa independent (black 
edges). 

GUILT BY ASSOCIATION 

Next we used the phosphatase interactome derived from the 
in vivo pull down experiment to ask whether the phosphatase 
interaction network could provide hints toward specific path- 
ways that are affected by phosphatase activity. To this aim, 
we performed a KEGG- pathway enrichment analysis of the 
ligands of each of the phosphatases, by using the Functional 
Annotation Tool, David (Huang da et al., 2009). The two 
phosphatases PTPN21 and PPP2CA and the PP2 scaffold sub- 
unit PPP2R1A were significantly associated to RTK signaling 
(Figure 3), in agreement with their established involvement in 
the modulation of EGF signaling by controlling the SRC and 
S6K kinases, respectively (Cardone et al, 2004; Carlucci et al, 
2010; Hahn et al., 2010). In addition our enrichment analy- 
sis reveals a statistically significant association of PPP3CA with 
cell differentiation signaling. This result is consistent with the 
report by Kao et al. that the differentiation of Schwann cells 
requires the activity of the PPP3CA phosphatase (Kao et al, 
2009). Similarly, we found that DUSP26 is significantly corre- 
lated with the DNA damage response. This conclusion is in 
accordance with the observations that DUSP26 inhibits the p53 
tumor suppressor function, by suppressing doxorubicin-induced 
apoptosis in human neuroblastoma cells (Shang et al., 2010). 
On the other hand, our experimental strategy led to the iden- 
tification of new biological processes that are controlled by 
these protein phosphatases (e.g., vesicular trafficking and cell 
metabolism). As shown in Figure 3, the interactors of both 
DUSP18 and PTPN3 were not significantly associated to any 
specific pathway. 

Next we looked for evidence that the proteins that copuri- 
fied with each phosphatase may form complexes. To this end, 
we queried the mentha database (Calderone et al., 2013), and 
looked for evidence of interactions between ligands of each bait 
phosphatase. The interactors of five of the eight phosphatases are 
linked by direct interactions. As illustrated in Figure 4A, we found 
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FIGURE 1 | Phosphatases affecting the RAS-PI3K pathway. (A) The results 
depicted in this figure are a subset of the results reported by Sacco et al. 
(2012a). The squares represents the main proteins that are involved in the 
PI3K and RTK/RAS signaling cascades. Green squares are the 
phosphorylated proteins that have been used as readouts in the screening. 
The phosphatase whose down regulation was shown to affect the readouts 
are represented as circles and are linked to readouts with red or green edges 
depending on whether their down-regulation negatively or positively affected 



the readouts. The intensity of the edge red or green color represents the 
trust in the reported interaction (number of concordant siRNA). The color of 
the phosphatase nodes is mapped to the phosphatase family according to 
the figure legend. (B) As previously described by Sacco et al., after 
phosphatase down-regulation the activation level of ERK and rpS6 was 
analyzed by automated immunofluorescence microscopy. Here, we report 
the ERK and rpS6 phosphorylation level after the down-regulation of some 
phosphatases hits. 



that DUSP26 copurifies with the serine/threonine kinase ATM, 
which, in response to genotoxic stress, phosphorylates the two 
Fanconi proteins FANCI and FANCD2, triggering the S-phase 
checkpoint activation (Taniguchi et al., 2002). A third DUSP26 
ligand, TEL02, which is a member, together with TTI1 and 



TTI2, of the TTT complex (Hurov et al., 2010) also interacts 
with ATM. 

PTPN21, on the other hand, interacts with the scaffold pro- 
tein GRB2, which associates with DNAJB11, DYNLL1, and UBR, 
suggesting that some of the identified interactors may copurify by 
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FIGURE 2 | Experimental strategy. Schematic overview of the experimental 
strategy applied to analyze the phosphatase interactome. Phosphatase 
transfected HeLa cells with or without TNFa were grown in media containing 
heavy isotope labeled lysine and arginine amino acids, while the control cells 
were grown in a medium containing natural amino acids. After lysis, 



phosphatases, and regulatory subunits were affinity purified by Strep 
purification and analyzed by mass spectrometry as described in Materials and 
Methods section. Phosphatase interactors were identified and quantified by 
their significant enrichment compared to the control and TNFa-induced 
alterations were quantified. 



indirect interactions (Figure 4C). Similarly we found that PTPN3 
interact with the mitochondrial ribosomal subunit ICT1 that 
binds GADD45GIP1 and POLRMT1 proteins (Figure 4B). 

As expected, the catalytic and scaffold subunit of PP2A share 
many interactors, confirming that such heterodimer forms differ- 
ent protein complexes that act on distinct substrates, by recruiting 
multiple regulatory subunits (Figure 4D). 

PTPN21 ASSOCIATES WITH THE SH3 DOMAIN OF GRB2 

Among all phosphatase-interaction partners, we focused on the 
newly discovered interaction between the scaffold protein GRB2 
and the tyrosine phosphatase PTPN21, both partners mapping 
to the RAS-PI3K signaling pathway. Cardone et al. reported that 



PTPN21 is recruited to mitochondria by binding the scaffold 
protein AKAP121 and that this interaction is essential for the 
phosphatase dependent dephosphorylation of the inhibitory 
tyrosine 527 of the SRC kinase (Cardone et al, 2004; Carlucci 
etal., 2010). 

GRB2 is an essential adapter protein consisting of two SH3 
domains flanking one central SH2 domain. The affinity purifi- 
cation assay results suggest that the GRB2-PTPN21 interaction 
is not likely to occur in a phosphorylation dependent manner, 
since it is not modulated by TNFa (Figure 3). However, phos- 
phoproteomics of both cancer and embryonic stem cells revealed 
that PTPN21 contains multiple tyrosine phosphorylated residues 
(Rikova et al, 2007; Guo et al, 2008; Brill et al, 2009). In order to 
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FIGURE 3 | The human phosphatase interactome. Phosphatases 
(squares) are linked to the experimentally identified interactors (circles) by 
edges. Edges are colored according to the functional relationships 
between the nodes they connect: interactions positively regulated by TNFa 
are in green; interactions negatively regulated by TNFa are red and TNFa 
independent interactions in black. Dashed lines represent interactions that 



have already been reported in the literature. Interactors are colored 
according to their functional association as revealed by our Kegg-Pathways 
enrichment analysis performed by the DAVID software. The phosphatase 
nodes are labeled according to the Kegg pathways that was significantly 
overrepresented in the phosphatase interactors and substrates 
(p-value < 0.005). 



map the GRB2-PTPN21 interaction to a specific GRB2 domain 
and assess whether such binding occurs in a phosphorylation 
dependent manner, the SH2 domain of GRB2 as well as its two 
SH3 domains were purified as GST fusion proteins and incubated 
with whole protein extracts co-transfected with Flag-PTPN21 
in presence or in absence of a constitutively active SRC kinase 
mutant (Y527F). As shown in Figure 5A, PTPN21 strongly binds 
the C-terminal SH3 domain of GRB2 and to a lesser extent the 
N-term SH3 domain, independently from SRC. On the other 
hand, the GRB2 SH2 domain does not interact with PTPN21. 
The analysis of PTPN2 1 protein sequence reveals that it contains 
a SH3 binding motif ( 564 RPPPPYPPPRP 574 ), whose sequence 
matches the GRB2 binding specificity described by Carducci et al. 
(2012). These results support the existence of a PTPN21-GRB2 
complex that is phosphorylation independent and likely occurs 
between the carboxy- terminal SH3 domain of GRB2 and PTPN2 1 
(Figure 5B). 



Next we asked whether the formation of the PTPN21 GRB2 
complex promotes the dephosphorylation of GRB2. For this 
purpose, GRB2 tyrosine phosphorylation was induced by trans- 
fecting HeLa cells with the constitutively active SRC kinase 
mutant (Y527F) in presence or in absence of Flag-PTPN21. As 
shown in Figure 5C, after cell lysis and immunoprecipitation 
with anti-GRB2, PTPN21 was found to associate with GRB2. 
However, when PTPN21 is overexpressed, GRB2 phosphoryla- 
tion, if anything, seems to be slightly increased, as revealed by 
probing the GRB2 protein with an anti-phospho tyrosine anti- 
body (Figure 5C). Thus, GRB2 is not a substrate of PTPN21 but 
may play a role in targeting PTPN21 to different substrates. 

A STRATEGY TO IDENTIFY NEW PHOSPHATASE SUBSTRATES IN 
GROWTH PATHWAYS 

Having obtained a high coverage interactome of the eight phos- 
phatases that affect the RAS-PI3K pathway we used it to develop 
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FIGURE 4 | Evidence for detection of multi-protein complexes. We 

queried the mentha protein interaction database to look for direct 
interactions among the phosphatase interactors. For five out of the eight 
baits we were able to retrieve information on direct interactions between 
the affinity purified interactors. Phosphatases are represented as yellow 
squares. DUSP26 (A), PTPN3 (B), PTPN21 (C), and PPP2CA (D) are linked 



to their substrates by different edges. As in Figure 1, continuous and 
dashed gray lines represent experimental and literature supported 
interactions, respectively. Black edges indicate direct interactions occurring 
among the phosphatase interactors. The color of the nodes is mapped 
according to the size of the putative complexes formed by the proteins 
according to the figure legend. 



a general strategy that could infer the direct target of these 
phosphatases. Phosphatase-substrate interaction is weak and 
transient, thus it is unlikely that substrates can be identified by co- 
immunoprecipitation. In fact none of the interactors identified 
in the affinity purification experiments are among the validated 
substrates annotated in the HUPHO and DEPOD databases (Li 
et al., 2013; Liberti et al., 2013). It has been reported that much of 
phosphatase substrate specificity, localization and activity is mod- 
ulated by the interaction with scaffold/regulatory proteins that 
target them to specific locations (Roy and Cyert, 2009; Sacco et al, 
2012b). We hypothesized that some of the interactors identified 
by our approach act as molecular bridges linking phosphatase to 
substrates participating in the RAS-PI3K pathway. For this reason, 



we made use of the PPI network downloaded from the mentha 
database (Calderone et al., 2013) to link phosphatase interactors 
to putative substrates in the RAS-PI3K pathway (Figure 6A). 
The strategy that we used is based on the following steps: 

(1) Draw a literature derived directed network of the RAS-PI3K 
pathway and identify the participating proteins as puta- 
tive targets of the "modulator phosphatases" (Supplementary 
Material, Table S4). 

(2) Combine the phosphatase interactors identified in the affin- 
ity purification MS experiment with the ones already 
described in the literature and reported in the mentha 
database (red and black edges, respectively, in Figure 6B). 
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FIGURE 5 | The PTPN21 phosphatase physically associates with the 
C-terminal SH3 domain of GRB2. (A) GST fusion proteins of the N-terminal, 
C-terminal SH3 domain as well as the SH2 domain of GRB2 were incubated 
with a protein extract of HeLa cells transiently transfected with Flag-PTPN21 
in presence or in absence of constitutively activated Y527F SRC kinase. 
Affinity-purified SH2 and SH3 domains ligands were separated by SDS-PAGE 
and transferred onto cellulose membranes. The blots were probed with 
anti-Flag (WB: a-4G10) and anti-phosphotyrosine (WB: a-4G10) antibodies. 
The cell lysate (input) and the sample affinity-purified with the GST protein 



(GST) were used as controls. (B) A schematic representation of the modular 
domain structure of GRB2 and PTPN21 proteins is represented. According to 
our data, the SH3 binding motif ( 564 RPPPPYPPPRP 574 ) of PTPN21 binds the 
GRB2 C-terminal SH3 domain. (C) HeLa cells were transiently co-transfected 
with Flag-PTPN21 and with constitutively activated Y527F SRC kinase 
expression plasmids. After cell lysis, whole protein extracts were 
immunoprecipitated with anti-GRB2 antibody. The membrane was probed 
with anti-GRB2 (WB: a-GRB2), anti-Flag (WB: a-Flag) and anti-phospho 
tyrosine (WB: a-4G10) antibodies. 



(3) Define paths in the protein interaction graph that connect 
each phosphatase to the proteins participating in a given 
pathway (here RAS-PI3K signaling). 

By this strategy, each interactor was linked to RAS-PI3K signaling 
proteins and a by a large number of possible paths. The result- 
ing complex graph was filtered according to the following rules 
(illustrated in Figure 6A): 

(1) Longer paths are filtered out. Only paths connecting tyro- 
sine phosphatases to protein members of the growth network 



with up to two "binding steps" are considered. For phos- 
phatases subunits that form holoenzymes with regulatory 
subunits such as PP2A and PPP3CA we allowed three binding 
steps. 

(2) The inferred substrate in the RAS-PI3K pathway has to be 
a kinase, a phosphatase or a scaffold protein whose dephos- 
phorylation controls either enzyme activity or the molecular 
association with other regulatory proteins. 

(3) The inferred substrate has to contain phosphorylation sites 
with defined functional roles as annotated in the Phosphosite 
database (Hornbeck et al, 2012). 
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FIGURE 6 | Inferring new phosphatase substrates. (A) Schematic 
representation of the multiple paths going from phosphatase to substrates. 
(B) The multiple paths going from phosphatases to substrates are 
represented as a graph. Nodes have different shapes according to their 
functional role: phosphatases are indicated as squares, bridge proteins as 
diamonds, modulators as hexagons and inferred substrate as circles. The red 
border outlines phosphatase substrates that have been already reported in 
literature. Solid black and red lines indicate physical interaction literature and 
experimentally supported, respectively, while black dashed line represent 



enzymatic interaction already described in literature. (C) HeLa cells were 
transiently co-transfected with expression plasmids expressing Flag-SHP2 
and a constitutively active mutant of the SRC kinase (Y527F) expression 
plasmids. After cell lysis, whole protein extracts were immunoprecipitated 
with anti-SHP2 antibody. The beads were washed with lysis buffer, and the 
immunoprecipitation (IP) was revealed with anti-SHP2 (WB: a-SHP2), 
anti-GRB2 (WB: ct-GRB2), anti-Flag (WB: a-Flag) and anti-phospho tyrosine 
(WB: a-4G10) antibodies. GRB2 which is an established ligand of SHP2 was 
used as a positive control. 
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(4) The phosphorylation sites of the inferred substrates have 
to be compatible with the nature of the phosphatase 
(Tyrosine phosphatase can only dephosphorylate tyrosine 
residues, ect). 

(5) The effect of the phosphatase induced dephosphorylation of 
the inferred substrate has to explain the phenotypic effect 
observed upon phosphatase down-regulation in the func- 
tional screenings (Sacco et al, 2012a). 

The result of this approach (Table S5) is illustrated in the filtered 
graph in Figure 6B. Remarkably, our strategy was validated by 
the recovery of phosphatase substrates already reported in the 
literature. For instance, Duan and Cobb already demonstrated 
that PPP3CA induces the MAPK activation by dephosphorylat- 
ing Thr401 in RAF1 (Duan and Cobb, 2010). In addition the 
inhibitory effect of PTPN3 on ERK phosphorylation was already 
reported by Han et al. (2000). Interestingly, both PPP3CA and 
PPP2CA phosphatases have been already described to be nega- 
tive modulators of autophagy (Magnaudeix et al., 2013; He et al., 
2014). Our approach enabled the identification of a new potential 
molecular mechanism that these two phosphatases may control to 
modulate autophagy. SIK3 and SQSTM proteins have been iden- 
tified by our affinity purification experiment as two novel inter- 
actors of PPP3CA and PPP2CA, respectively. In our approach, we 
propose that SIK3 and SQSTM proteins act as bridge to connect 
PPP3CA and PPP2CA phosphatases to the autophagy marker 
MLP3A (LC3A). This observation suggests that our approach can 
be used to propose new potential molecular mechanisms link- 
ing a phosphatase to an established biological process. This graph 
links phosphatase to putative adapter and to putative substrates. 
In principle depending on the available information one can use 
it (1) to infer new substrates starting from a consolidated PPI 
or (2) to validate molecular bridges that target a phosphatase to 
an established substrate. In the two following paragraphs we will 
demonstrate these strategies in two specific cases. 

SHP2 CAN BE DEPH0SPH0RYLATED BY DUSP18 

DUSP18 was shown by our screening to negatively regulate the 
RAS pathway. The graph in Figure 6B indicates that the regula- 
tory protein that is closest to DUSP 1 8 in the RAS pathway is SHP2 
and that DUSP 18 and SHP2 are connected by catalase. Indeed 
it has been shown that the SH2 domains of SHP2 bind tyro- 
sine phosphorylated catalase (Yano et al., 2004), and catalase was 
recovered as a DUSP 18 interactor in our approach. We can there- 
fore picture catalase acting as a bridge linking the phosphatase 
to its putative target. To test this hypothesis, HeLa cells were 
transiently co-transfected with Y527F SRC kinase, to enhance 
phosphorylation, and Flag-DUSP18. As shown in Figure 6C, after 
cell lysis and endogenous immunoprecipitation with anti-SHP2, 
DUSP 18 was found to associate with SHP2 only in SRC trans- 
fected cells. This data is compatible with the model whereby 
the SH2 domains of SHP2 bind tyrosine phosphorylated catalase 
which in turn binds to DUSP 18. In addition the over-expression 
of DUSP 18 induces SHP2 dephosphorylation, without affect- 
ing its association with GRB2. Since it has been shown that the 
C-terminal tyrosine residues of SHP2 bind GRB2, this result 
suggests that DUSP 18 likely dephosphorylates the Tyr62 and 



Tyr63 residues. Although the biological relevance of the inferred 
dephosphorylation needs to be proven in more physiological 
conditions, this result shows that DUSP 18 has the potential to 
dephosphorylate SHP2 as inferred by our approach. 

SCRIB ACTS AS A BRIDGE TO TARGET DUSP26 TO ERK 

Knock down of DUSP26 by siRNA negatively affects the activa- 
tion of ERK (Sacco et al., 2012a). This is in agreement with the 
ability of DUSP26 to inhibit cell proliferation in epithelial cell 
lines (Hu and Mivechi, 2006; Patterson et al, 2010). Consistent 
with a role as tumor suppressor, DUSP26 is down-regulated, 
in several human cancer cell lines, as well as in some primary 
tumors (Tanuma et al, 2009; Patterson et al., 2010). However, 
DUSP26 is not able to directly bind ERK to dephosphorylate it 
(Hu and Mivechi, 2006; Patterson et al., 2010) suggesting the exis- 
tence of a molecular bridge The heat shock transcription factor 
Hsf4b, a substrate of ERK, was proposed as a possible bridge to 
link DUSP26 to ERK (Hu and Mivechi, 2006). Similarly, more 
recently, the adenylate kinase 2 was proposed to be a bridge that 
directs DUSP26 to dephosphorylate FADD (Kim et al, 2014). 

Our approach identified SCRIB as a potential bridge that 
would modulate the de-phosphorylation of ERK by DUSP26. 
SCRIB is an adapter protein that was recently suggested to down- 
regulate ERK by binding and activating the phosphatase PPI 
gamma (Nagasaka et al., 2013). We propose here that SCRIB 
may also promote the de-phosphorylation of ERK by DUSP26. 
SCRIB directly binds to ERK through two KIM motifs and reg- 
ulates its activation and nuclear translocation (Nagasaka et al., 
2010). The protein contains four PDZ domains (Figure 6A). The 
C-terminus region of DUSP26 contains an atypical motif for PDZ 
binding L-D/E-<t>, where <t> is a hydrophobic residue (Tonikian 
et al., 2008). Thus, we asked whether the binding of SCRIB to 
DUSP26, as identified in our affinity purification experiment, 
could be mediated by any of the SCRIB PDZ domains. To this 
end we performed a GST pull down experiments by affinity puri- 
fying extracts of E. colt cells expressing HIS-tagged DUSP26 with 
GST fusion of SCRIB PDZ domains (Figure 7A). Only the fourth 
PDZ domain of SCRIB was able to bind DUSP26. The binding 
was confirmed by co-immunoprecipitation assay, after cotrans- 
fecting HA-SCRIB and Flag-DUSP26 in H1299 cells. As shown in 
Figure 7B, co-immunoprecipitated SCRIB was detected by west- 
ern blotting with an anti-SCRIB antibody. Similarly SCRIB was 
immunoprecipitated with anti-HA antibody and the presence of 
DUSP26 was revealed by western blotting with anti-DUSP26 anti- 
body Figure 7C. These data suggest that SCRIB could direct the 
phosphatase activity of DUSP26 toward ERK as suggested in the 
cartoon in Figure 7D. 

DISCUSSION 

Although protein phosphorylation has been considered as a 
key post-translational mechanism controlling a variety of phys- 
iological processes and a number of reports have contributed 
to describe the phosphatase interaction network, a compre- 
hensive characterization of phosphatase substrates is still miss- 
ing (Goudreault et al, 2009; Breitkreutz et al, 2010; Skarra 
et al., 2011; Couzens et al., 2013). Recently we have reported 
an unbiased siRNA screening aimed at identifying phosphatases 
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FIGURE 7 | DUSP 26 interacts with the adapter protein SCRIB. (A) GST 

Pull Down of full-length DUSP26 by GST constructs fused to PDZ 1-2, PDZ3, 
PDZ4, or all the four PDZ domains (aa 728-1630) of SCRIB fused to GST 
Black Lines under the schematic representation of the domain structure of 
SCRIB indicate the protein regions that were fused to GST. Two controls 
were added: GST alone and a PDZ containing protein of similar length 



(PDZK1, lane "control"). (B,C) Co-immunoprecipitation of SCRIB and 
DUSP26. H1299 cells were transfected with HA-SCRIB, Flag-DUSP26 or both 
plasmids as indicated. Lysates were immunoprecipitated with anti-Flag (B) or 
with anti-HA (C) and detection was performed with anti-DUSP26 or 
anti-SCRIB antibodies. (D) Cartoon picturing the proposed role of SCRIB as 
adapter protein to mediate DUSP26 dephosphorylation specificity. 
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controlling key growth pathways in cancer cells (Sacco et al., 
2012a). Combining the siRNA screening results with modeling 
techniques, we were able to map phosphatases on specific nodes 
of the growth signaling model. However, our approach only iden- 
tified phosphatases modulating the growth pathway but did not 
enable us to link phosphatases to specific substrates. 

For this purpose, we set up to develop an experimental strat- 
egy that combines the functional information obtained with the 
siRNA screening and PPI network context information. We first 
enriched the literature derived interactome of six phosphatases 
and two phosphatase accessory subunits by affinity purification 
experiments of phosphatase complex followed by quantitative 
mass spectrometry based proteomics in cancer cells stimulated 
with TNFa. By this approach we were able to recapitulate most of 
the interactions occurring between the catalytic, scaffold, and reg- 
ulatory subunit of the PP2A holoenzyme, confirming the robust- 
ness of our approach. The resulting interactome is completely 
connected, since each phosphatase shares at least one ligand 
with one of the remaining phosphatases. For instance the tyro- 
sine phosphatase PTPN21 and the catalytic subunit of the serine 
threonine phosphatase PP2A share a common group of interac- 
tors, mainly involved in controlling cell metabolism. We observed 
that the phosphatase interactome is largely insensitive to stimula- 
tion with TNFa, suggesting that these interactions may be either 
constitutive or triggered by other types of stimuli. For instance, 
while the DUSP26-ATM interaction is not modulated by TNFa, 
we show that nutrients and amino acids deprivation increases 
the binding (Supplementary Material, Figure S1A), suggesting 
that these proteins may play a role in controlling the autophagy 
process. Indeed, we have previously shown that the siRNA inter- 
ference of DUSP26 results in a decrease of the autophagy marker 
LC3, while much evidence suggest that the ATM kinase pro- 
motes the autophagy induced by ionizing radiation and ROS 
(Liang et al, 2013; Tripathi et al., 2013). In addition, as shown 
in Figure 3, about 50% of the PPP3CA interactions are nega- 
tively modulated by TNFa, including the binding to its activator 
subunit calmodulin. This result suggests that the TNFa stimula- 
tion may have an inhibitory role on PPP3CA activity. However, 
Fernandez et al. have recently demonstrated that in reactive, 
but not in quiescent astrocytes, PPP3CA dephosphorylates the 
transcription factor Foxo3 in response to TNFa, suggesting that 
depending on the cells type, this phosphatase may have opposite 
functions (Fernandez et al., 2012). 

Interestingly our experimental approach enabled us to iden- 
tify a novel interaction between the scaffold protein GRB2 
and the tyrosine phosphatase PTPN21. Here, we report that 
PTPN21 binds the C-terminal SH3 domain of GRB2 in vitro, 
but it does not dephosphorylate its phosphotyrosine residues. 
Indeed our affinity purification experiment failed to identify 
known phosphatase substrates that had already been described 
in the literature. This observation is not surprising if we 
consider that phosphatases rapidly dephosphorylate the sub- 
strate and the phosphatase-substrate interaction is so tran- 
sient and weak that coimmunoprecipitation-based approaches 
likely fail to identify phosphatase substrates. In addition, while 
most protein kinases recognize a specific amino acid motifs 
on their targets, phosphatase substrates specificity is weaker 



and mainly based on the interaction with regulatory subunits 
(Roy and Cyert, 2009). 

To infer new phosphatase substrates, we have here outlined 
a combined experimental-bioinformatic strategy based on the 
integration of the phosphatase interactome with network context 
information, extracted from the mentha PPI database (Calderone 
et al, 2013). Although this approach lead us to recover some of 
the phosphatase-substrate relationships already described in lit- 
erature, we are aware of some relations that are missed by our 
approach [e.g., the RAF1 dephosphorylation by PPP2CA (Dent 
et al., 1995)]. These failures can be explained by several factors: 
(1) some interactions may be cell type dependent or rely on spe- 
cific stimulations; (2) some phosphatase partners may have very 
low level of expression that remains undetected in our affinity 
purification experiments and (3) some PPI relations may have not 
been reported yet or may have not been annotated in mentha. In 
addition we want to stress that we used rather stringent filtering 
criteria to reduce the total number of inferred phosphatase- 
scaffold-substrate complexes. This might increase the chance of 
missing already validated enzyme-substrate relationships or of 
identifying new interesting regulation mechanisms. If desirable, 
these criteria can be relaxed at the cost of increasing the noise of 
false positives. 

In essence our method combines functional information with 
the interactome and analyses the resulting graph to identify paths 
between phosphatases and putative substrates. By this approach 
new substrates may be inferred or alternatively proteins that 
form molecular bridges between the phosphatase and the sub- 
strates can be identified. To assess the robustness and reliabil- 
ity of our strategy, two specific cases were analyzed. Firstly we 
demonstrated that DUSP18 induces SHP2 dephosphorylation. 
Our siRNA screening revealed that DUSP18 negatively controls 
ERK phosphorylation (Sacco et al., 2012a). This is consistent 
with SHP2 being a positive modulator of the MAPK signaling 
(Cai et al, 2002). Here we infer that catalase acts as a bridge 
to enable the DUSP18 mediated de-phosphorylation of SHP2. 
Although our approach does not identify the specific SHP2 tyro- 
sine residues dephosphorylated by DUSP18, we demonstrated 
that the C-terminal residues involved in the GRB2 interaction are 
not targeted by the phosphatase (Figure 5). DUSP18 may neg- 
atively controls the MAPK signaling by dephosphorylating and 
inactivating SHP2. Finally, we demonstrated that SCRIB acts as 
a bridge to mediate the dephosphorylation of ERK by DUSP26 
(Figure 7B). DUSP26 is a poorly characterized dual specificity 
phosphatase whose negative regulation of the MAPK signaling 
has been already reported. 

Taken together these observations show that the combina- 
tion of the topological information contained in the phosphatase 
interactome with functional information obtained by siRNA 
screening can be valid tool to infer new phosphatase substrates 
and modes of targeting. 

MATERIALS AND METHODS 
ANTIBODIES AND REAGENTS 

Anti-HA, anti-FLAG and anti-Flag Ml agarose beads and anti 
DUSP26 were from Sigma; anti-SHP2 and anti SCRIB were 
from Santa Cruz Biotechnology; anti-GRB2 and anti-4G10 was 
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from Upstate Biotechnology, Inc. Peroxidase-conjugated anti- 
rabbit, anti-mouse and anti-goat secondary antibodies were 
from Jackson ImmunoResearch. PPP2CA, PTPN3, PTPN21, 
DUSP26, PPP2R3C encoding plasmids were purchased from 
OpenBiosystem. DUSP18, PPP2R1A, and PPP3CA constructs 
were kindly provided by Marc Vidal. Phosphatase cDNAs were 
cloned in pDNOR vector (Invitrogen) and cloned in the SF-TAP 
plasmid by using the Gateway Recominant Cloning Technology 
from Invitrogen. The cDNA of DUSP26 was also cloned in Pet28 
and PC-DNA plasmids. HA-DLC1 was kindly provided by Prof. 
Cecconi. The cDNA encoding SRC Y527F was cloned in pSGT 
(Gonfloni et al, 1997). HA-SCRIB, PDZ3-GST, and PDZ4 -GST 
were a generous gift of L. Banks. Construct containing human 
SCRIB PDZ1-2 and 1-4 (aa 728-1630) were cloned in pGex2TK. 

CELL CULTURE 

Cells were maintained in a humidified atmosphere at 37°C and 
5% C02 in Dulbecco's modified Eagle's medium (Invitrogen), 
supplemented with 10% fetal bovine serum (Sigma) and 0.1% 
penicillin/streptomycin (Invitrogen). For SILAC experiments, 
SILAC DMEM (PAA, Pasching, Austria) deficient of L-Lysine and 
L-Arginine, supplemented with 10% (v/v) dialyzed fetal bovine 
serum (FBS; PAA, Pasching, Austria), 50 units/ml Penicillin, 
0.05 mg/ml Streptomycin and 0.55 mM lysine, 0.4 mM arginine 
was used. Light labeled medium was supplemented with u Ce, 
14 N2 lysine and 12 Cs, 14 N4 arginine, medium labeled medium 
with 4.4.5.5-D4-L-Lysine and 13 C6- 14 N4-L-Arginine and heavy 
labeled medium with 13 C 6 15 N 2 -L-Lysine and 13 C 6 15 N 4 -L- 
Arginine. Proline was added to a final concentration of 0.5 mM 
to prevent arginine to proline conversion (Bendall et al, 2008), 
which could impair the quantification. All amino acids were pur- 
chased from Silantes. Human epithelial carcinoma (HeLa) cells 
were purchased from the ATCC. HeLa cells were transfected with 
Lipofectamine 2000 (Invitrogen) according to manufacturer's 
protocol. 

AFFINITY PURIFICATION OF PROTEIN COMPLEXES 

For one step Strep purifications, SF-TAP tagged proteins 
and associated protein complexes were purified essentially as 
described earlier (Gloeckner et al., 2007; Boldt et al., 2011). 
HeLa cells, transiently expressing the SF-TAP tagged constructs 
or SF-TAP alone as control were either stimulated with 50 ng/ml 
TNFoi or mock treated. They were next lysed in lysis buffer (con- 
taining 150 mM NaCl, 50 mM Tris-HCl, 1% Nonidet P-40, and 
0.25% sodium deoxycholate, protease inhibitor cocktail (Roche) 
and phosphatase inhibitor cocktails II and III (Sigma- Aldrich), 
for 20min at 4°C. After sedimentation of nuclei at 10,000 x g 
for lOmin, the protein concentration of the lysates were deter- 
mined by a Bradford assay before equal amounts of the cleared 
lysates were transferred to Strep-Tactin-Superflow beads (IBA) 
and incubated for 1 h before the resin was washed three times with 
wash buffer (TBS containing 0.1% NP-40, phosphatase inhibitor 
cocktail I and II). The protein complexes were eluted by incuba- 
tion for lOmin in Strep-elution buffer (IBA). Following elution, 
the corresponding samples were combined. The combined sam- 
ples were concentrated using lOkDa cut-off VivaSpin 500 cen- 
trifugal devices (Sartorius Stedim Biotech) and pre-fractionated 



using SDS-Page and in-gel tryptic cleavage as described elsewhere 
(Gloeckner et al, 2009). 

MASS SPECTROMETRY AND DATA ANALYSIS 

LC-MS/MS analysis was performed on an Ultimate3000 nano 
RSLC system (Thermo Fisher Scientific) coupled to a LTQ 
Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) by 
a nano spray ion source. Tryptic peptide mixtures were automat- 
ically injected and separated by a linear gradient from 5 to 40% 
of buffer B in buffer A (2% acetonitrile, 0.1% formic acid in 
HPLC grade water) in buffer A (0.1% formic acid in HPLC grade 
water) at a flow rate of 300 nl/min over 90 min. Remaining pep- 
tides were eluted by a short gradient from 40 to 100% buffer B 
in 5 min. The eluted peptides were analyzed by the LTQ Orbitrap 
Velos mass spectrometer. From the high resolution MS pre-scan 
with a mass range of 300-1500, the 10 most intense peptide ions 
were selected for fragment analysis in the linear ion trap if they 
exceeded an intensity of at least 500 counts and if they were at least 
doubly charged. The normalized collision energy for CID was set 
to a value of 35 and the resulting fragments were detected with 
normal resolution in the linear ion trap. The lock mass option 
was activated, the background signal with a mass of 445.12002 
was used as lock mass (Olsen et al., 2005). Every ion selected 
for fragmentation, was excluded for 20 s by dynamic exclu- 
sion. For SILAC experiments, all acquired spectra were processed 
and analyzed using the MaxQuant software (Cox and Mann, 
2008) (version 1.0.13.13) and the human specific IPI database 
version 3.52 (http://www.maxquant.org/) in combination with 
Mascot (Matrix Science, version 2.2). Cysteine carbamidomethy- 
lation was selected as fixed modification, methionine oxidation 
and protein acetylation were allowed as variable modifications. 
The peptide and protein false discovery rates were set to 1%. 
Contaminants like keratins were removed. Proteins, identified 
and quantified by at least two unique peptides were considered 
for further analysis. The significance values were determined by 
Perseus tool using significance B. Those proteins whose ratio 
was greater than 1.9 and significance B was lesser than 0.1 were 
considered significantly enriched. 

PULL-DOWN ASSAY 

After 24 h of transfection, confluent HeLa cells were washed with 
ice-cold PBS and lysed in RIPA buffer (150 mm NaCl, 50 mm 
Tris-HCl, 1% Nonidet P-40, 0.25% sodium deoxycholate) sup- 
plemented with 1 mm pervanadate, 1 mm NaF, protease inhibitor 
mixture 200 x (Sigma), inhibitor phosphatase mixture I and II 
lOOx (Sigma). The samples were kept on ice for 30 min and 
centrifuged at 15,000 rpm at 4°C for 30 min. The supernatant 
was collected, and the total amount of protein was determined 
by Bradford colorimetric assay (Bio-Rad). The whole cell lysates 
were incubated with 50 |ig of the indicated GST fusion protein at 
4°C for 1 h. Thus, glutathione-Sepharose 4B beads were blocked 
by incubating with 3% bovine serum albumin with rocking at 
4°C for 1 h, and then after centrifugation for 3 min at 4000 x 
g, at 4°C, the dry beads were bound to lysates mixed with GST 
fusion proteins at 4°C for 1 h. The supernatant was discarded by 
centrifugation, and the beads were washed six times with lysis 
buffer for 3 min at 4000 x g, at 4°C, and then the dry beads 
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were resuspended in SDS sample buffer, boiled and analyzed by 
SDS-PAGE and Western blotting on nitrocellulose membrane. 

IMMUNOPRECIPITATION AND IMMUNOBLOT ANALYSIS 

HeLa cells were lysed as described previously. The whole cell 
lysates were incubated with anti-Flag antibody conjugated to 
Sepharose beads over-night at 4°C. The beads were washed with 
lysis buffer, and the immunoprecipitated proteins were sepa- 
rated by SDS-PAGE, transferred onto a nitrocellulose membrane, 
and immunoblotted with antibodies. The immunoreactions were 
visualized using ECL detection system (Amersham Biosciences). 
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SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found online 
at: http://www.frontiersin.org/journal/10.3389/fgene.2014. 
00115/abstract 

Figure S1 | Validation of some of the newly identified phosphatase 
interactions (A) HeLa cells were transiently transfected with Flag-DUSP26 
expression plasmid. Twenty-four hours post transfection, cells were 
serum and amino acids starved for 1 h or left untreated and then lysed. 
Whole protein extracts were immunoprecipitated with anti-Flag antibody 
to purify the DUSP26 phosphatase. The membranes were probed with 
anti-ATM (WB: a-ATM) and anti-Flag (WB: a-Flag) antibodies. (B) HeLa 
cells were transiently co-transfected with Flag-PTPN21 and with HA-DLC1 
expression plasmids. Twenty-four hours post transfection, cells were 
serum and amino acids starved for 1 h or left untreated and then lysed. 
Whole protein extracts were immunoprecipitated with anti-Flag antibody 
to purify PTPN21 phosphatase. The membranes were probed with 
anti-HA (WB: ot-HA), anti-Flag (WB: a-Flag) and anti-GRB2 antibodies. (C) 
HeLa cells were transiently co-transfected with Flag-PPP2CA and with 
HA-DLC1 expression plasmids. Twenty-four hours post transfection, cells 
were serum and amino acids starved for 1 h or left untreated and then 
lysed. Whole protein extracts were immunoprecipitated with anti-Flag 
antibody to purify PTPN21 phosphatase. The membranes were probed 
with anti-HA (WB: a-HA) and anti-Flag (WB: a-Flag) antibodies. 

Table S1 | Protein groups identified by mass spectrometry based 
proteomics of phosphatase pull-down are reported with protein 
quantification, number of peptides and intensities. 

Table S2 | After statistical analysis, for each phosphatase the 
corresponding interactor is reported. In the "Phosph-SF" column, the 
intensity value of each interactor in phosphatase transfected cells was 
divided by its intensity in not transfected cells (Control). In the 
"PhosphTNF-Phosph" column, the intensity value of each interactor in 
cells over-expressing the phosphatase and stimulated with TNFa was 



divided by its intensity in transfected unstimulated cells. Finally in 
"PhosphTNF-SF" column, the intensity value of each interactor in cells 
over-expressing the phosphatase and stimulated with TNFa was divided 
by its intensity in not transfected cells. For each ratio, the corresponding 
significance B is reported. 

Table S3 | List of common contaminants was collected from the literature. 
Table S4 | Experimental data describing the functional relationships 
between signaling proteins in the pathways of interest were collected 
from the literature (PMID column). Each enzyme-substrate relationship is 
described as activating (1) or inhibitory (-1). For each protein, Uniprot ID 
and gene name have been reported. 

Table S5 | Experimental and literature extracted binary interactions, 
describing the paths from a phosphatase to its target in the RAS-PI3K 
pathway. 
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